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Abstract. In this note a general result is proved that can be used to evaluate exactly 
a class of highly oscillatory integrals. 



1. Introduction 

The first of the ten $100,100-digit challenges [H [5], proposed to calculate the inte- 
gral Jq x^^ cos{x^^ log x) dx to ten digits. This was a real challenge since the integrand 
oscillates infinitely often inside the interval of integration. 

In this note, a general result is proved that will allow us to determine exactly the 
value of some highly oscillating integrals. To give you the fiavour of what we will prove, 
here is one of our results : 

de _ ( + 3 - 2y/2 

Jo l + sin2(tan^) ~ 2^2 \^e2 -3 + 2^2 
where the graph of the integrand is depicted in the following Figure [TJ 




Figure 1. The graph of the function x ^ , , . — r on TO, f). 

° l+sin (tanx-j L '2/ 

In section [2] we will prove our main results and in section [3] we will give some detailed 
examples and applications. 
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2. The Main Results 

First, let us set the framework of our investigation. Our starting point will be a 
27r-periodic locally integrable function / G L^(T). For a nonzero real number A, we are 
interested in the evaluation of the integrals: 

Mf)= T^de=— /e(Atanx)dx, (2.1) 



A2 + ^2 U| 



and 



J7a(/) = ^ ^^^de = J^^ fo{\X\tanx)tanxdx, (2.2) 

where /e and fo are, respectively, the even and the odd components of /. They are defined 
on M by the formulae: 

m^imitu fM^rn^. (2,3) 

On the other hand, the Cauchy principal value C\{f) of the integral dt, plays 
an important role in this study. Recall that C\{f) is defined as follows: 

A(/)= lim r/^dt (2.4) 



In the next proposition, we prove the convergence of the integrals defined by fl2.ll) 
and (12. 2p . and we describe their relationship to the principal value defined by (12.41) . 

Proposition 2.1. Consider f G L^(T), and a nonzero real A. 
f°° fit) 

i. The integral / — — ■ dt is absolutely convergent. 

Jo t + X 

f°° tfit) 

ii. If Jjf{t)dt = 0, then the integral / ^2 _[_ ^2 ^"^ convergent. 
Hi. The principal value C\{f) is well defined and 

£,(/) = 2 (J,(/)- zAX,(/)). (2.5) 
where Xx{f) and J\{f) are defined by (12. ip and (12. 2p . respectively. 
Proof. Indeed, the first point is easy since, for /c > 0, we have 

''''''^\mLdt= r dt<^r\f{t)\dt 

So 



and dl]) follows. 
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To prove dH]) we consider the function F defined by F{x) = f{t) dt. Since / 
belongs to L^(T) and jjf{t)dt = we conclude immediately that F is continuous and 
27r-periodic. In particular, F is bounded on M. 

An integration by parts shows that, for x > 0, we have 

,f{t)dt = ——F{x)- -——-F{t)dt 



Indeed, this version of "integration by parts" is a dirct application of Fubini's theorem, 
see for instance [3l Theorem 5.2.3]. Now, the following inequality 



Vt > 0, 



M , , 

- A^rp^ where M = supjj 



(A2 + t2)i 

jj^-^^^yF{t)dt. This implies the 

f°° t 

convergence of / — r- /'(t)(it and shows that 

Jo \^ + t^^' 



. POO .2 _ \2 



On the other hand, for a > 0, we have 



a 



t + l\ J_„ A2 + t2 



because t tfeit) — iXfoit) is odd, and t h-)- tfo{t) — iX fe(t) is even. 

Finally, noting that fj, fo(t)dt = since is odd, and using @ and dH]), we obtain 
(12. 5p by letting a tend to +oo, and the Proposition follows. □ 

Using the notation of the preceding proposition, and the parity of A T\{f) and 
A I—)- J7a(/), we obtain the following corollary that reduces the determination of I\{f) and 
Jxif) to that of Cxif). 

Corollary 2.2. Using the notation of Proposition [K7\ we have 

Uf) = {^{^x{f)-^-x{f)), (2.6) 
Jx{f) = \{C,{f) + C^,{f)). (2.7) 

Thus, the question is reduced to determining Cx{f)- The next proposition expresses 
differently this quantity. 
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Proposition 2.3. For f G L^(T), and a nonzero real X, we have 



1 Z-^" ft + iX 







£,(/) = - / cot ( ] f{t)dt (2.8) 



where C\{f) is defined by (12.41) . 

Proof. Recall that, (see [21 pages 187-190], 



n-l 



cot 2; = - + > + = hm > 



z ^ \ z — nn z + Tin I n^oo ^ z + Tik 

n=l ^ ' k=—n 

with normal convergence on every compact set K contained in C \ vrZ. Applying this to 
the compact segment K = : t G [0,27r]}, and recalling that / G L^iT), we conclude 

that 

Jo \ 2 J^^' ™^t^„V-^o t + 2nk + iX 



t + tX 



= 2 lim V f / 

— ,t^„ 

= 2 lim I^dt, 
= 2£a(/). 

Which is the desired conclusion. □ 

The following lemma, will be of interest in formulating our main theorem. 

Proposition 2.4. For f G L^{T), and a nonzero real X, we have 

A(/) = m5^(sgn(n) - sgn(A))e-l"^IC„(/), (2.9) 



where sgn(x) is the sign of x if x 7^ and sgn{0) = 0, and {Cnif))^^^ are the exponential 
Fourier coefficients of f . 

Proof. For a nonzero real A, consider the continuous, 27r-periodic, function g\ defined by 

g,it) = cot (^-^] . (2.10) 



2 

The exponential Fourier coefficients {Cn{gx))nez of gx are given by 

CM = ^(sgn(A) - sgn(n))e-l"^l (2.11) 



Indeed, let us consider two cases 
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For A < 0, we have 



= 1 + J22e 

n=l 

And, for A > 0, we have 

, , e-^-** + 1 



Q—X—it ]^ ]^ g— A— ii 

oo —1 



X^int 



n=l n=— oo 

Combining these two points proves (12. lip . 

Now, since |C„(/)| < ||/||-^ for every n G Z, we see that 

J2\Cn{9x)\\Cn{f)\ <+00. 

Thus, the continuous, 27r-periodic, function gx* f is equal to its Fourier series expansion. 
In particular. 



9x * /(O) = ^ r f{t)9x{-t) dt = J2 Cn{9x)CnU) 
That is, by (ETTj) . 

^ r fit) cot (^-^) rft = -^J]((sgn(A)-sgnH)e-l"^l)C7„(/). 

Finally, using Proposition 12. 3[ we obtain f l2.9p . □ 

Now, combining Corollary 12. 2[ Proposition 12.11 and Proposition 12. 4[ we come to our 
main Theorem. 

Theorem 2.5. Consider f G L^(T) and a positive real X, then the following two integrals 
are convergent 

at ana / at 



A2 + t2 X2 + t2 

where fe and fo are, respectively, the even and odd components of f defined by (\2.3\) . 
Moreover, 

r^rft = |^5^sgn(n)e-H^C„(/). (2.13) 
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Corollary 2.6. Consider f G L'^{T) and a positive real X. 

:) = 

2 



i. If f is even, and its Fourier series expansion is S[f]{x) = — + ^^a„cosrix, then 



n>l 



oo 







A2 + t2 2A \ 2 



n=l 



n>l 



a. If f is odd, and its Fourier series expansion is S[f]{x) = ^^6„smna;, then 

r^dt=^Tbe-- 



There is a case where Theorem [23] takes a more practical form. So, let us change our 
point of view. Let G : Q — > C be an analytic function on a domain fl, that contains the 
closed unit disk D{0, 1). For a positive real number A, we are interested in determining 
the value of the integrals : 

9cio) _ 1 r/' 



.^U_d9=-j^ g,{Xtanx)dx, (2.14) 



and 



Ogsie) 



tt/2 



Jx{G) = / -— — 7^ d6 = / gs{Xtanx)tanxdx, (2.15) 
Jo X^ + 9^ Jq 

where the real variable functions Qc and Qs are, respectively, the "cos" and the "sin" 
components of the function 6 i— > G(e*^), defined on M by the formulae: 

G(e'^) + G(e-^) G(e-^) - G(e-^) 
^c(^) = and gs{6) = . (2.16) 

The next theorem gives the answer to this question : 

Theorem 2.7. Let G be an anlytic function on a domain Q that contains the closed unit 
disk, then, for every positive real X, one has 



J^,.^l^Gie-% (2.17) 



and 



A2 + x2 2 
where Qc and Qs are defined by fl2.16p . 



" dx = ^ (Gie-^) - G(0)) , (2.18) 
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Proof. Indeed, by assumption, there is a power series Yl'^=o ^n-^" with radius of conver- 
gence p > 1 such that 

oo 

VzgC, \z\<p =^ G(2) = J]a„z" (2.19) 

n=0 

So, we can define a continuous, 27r-periodic function / by f{d) = G(e*^). Comparing 
(12.1 6p and (12. 3p we see that fe{d) = gdO) and fo{0) = Wsi^)- On the other hand, from 
(I2.19P we conclude that the Fourier coefficients of / are given by C„(/) = a„ if n > and 
Cn{f) = if n < 0. Thus, Theorem 12.71 proves 

n>0 

"^0 n>0 

Which is the desired conclusion. □ 



3. Examples and Applications 

Example 3.1. 

Let X be a positive number, and consider the analytic function G defined on the 
domaine Q = C \ {e^^}, by 

p2x _|_ ^ 

= (3-1) 

Here, it is straightforward to check that 

Qf^^ie^ ^ sinh2x + ising 2) 
cosh 2x — cos 6 ' 



so 



cosh 2x — cos 9 
Using Theorem 12.71 we see that, for A > 0, we have 



sinh 2x 

'^^^^ = cosh 2x- cos ^' ^'-'^ 
9sie) - ,3^.?^. (3.4) 



A2 + ^2 2A e2^'-e 
Thus, from (13. 3p . and after making the change of variables 6 26 and A ^ 2A, we get 

1 de TT e2^ + e-2^ 



/o cosh 2x- cos 20 A2 + 02 2Asinh2x e^^ - e-^^ 
or, equivalently 

1 de _ TT e^^ + e-^^ 

sinh^ X + sin^ 6* A^ + 6*2 ~ A sinh 2x e2^ - e-2^ 



(3.6) 
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This can be expressed as follows : For positive A and fi, we have 

(3.7) 



In particular, choosing A = /i = 1 we find the integral (II .ip that we used to introduce our 
discussion. 

On the other hand, using Qs from fl3.4l) . and the second part of Theorem 12. 7[ we find 
that, for A > 0, we have 

^S-.^. (3.8) 



/o (cosh 2x - cos 6) (A^ + 6^) e2^+^ - 1 
or, equivalent ly, for /i > 1 and A > : 



d9 = , . (3.9) 



(/i-cose)(A2 + ^^2) eA(^+y;;^^)_i- 

The integrals in this example are to be compared with [3 formulae 3.792(10) and 3.792(13), 
page 450]. 

Example 3.2. 

In our second example we consider the even function / G L^(T) defined by 

/(t) = -ln|sin(t/2)| = -^lnsin2(t/2), for t ^ 27rZ. (3.10) 

It is well-known, (see for instance [lOi, Chapter 3, $ 14],) that / has the following Fourier 
series expansion 



Sm) =\n2 + J2 



CO 

COS nt 



n 

n=l 



Thus, applying Corollary 12.61 we obtain that for every A > we have 

/■~ln(sin2(t/2)) , TT / ^e-"^\ vr „ 
L ^A^^^--A(^^2 + X:^)^--(ln2-ln(l-e-)). 



The change of variables t ^ 2t and A ^ 2A yields the following result: 
, , r^ln(smH) , vr, fl-e''^\ 

Our next example is a generalization of an old problem. 
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Example 3.3. A Generalization of A Problem of Narayana Aiyar. 

Let ti,t2, . . . ,tn and a be positive real numbers, such that 

< ti < t2 < . . . < tn < a. 
Consider the meromorphic function G defined by 

G(z) = ; (3.12) 

[a — tiz){a — t2z) ■ ■ ■ [a — tnZ) 

Clearly, G is analytic in the domain r2 = C\ |^:l<A;<n| that contains the closed 
unit disk. 

For a given real 9, let 

0fc(6') = arctan ( ^^^^^^ ] ^ p^{Q^ = .ja^ _ 2tfcacos6' + t^. (3.13) 
\a — tk cos o J V 

To simplify the notation, we will simply write (pk and pk to denote (j)k{9) and pk{6) 
respectively. It is clear that 

a - tke'*^ = pke-"^\ for 1 < A; < n. (3.14) 

Thus, 



e 



i(4>l-\ h0n) 



G{e'') = . (3.15) 

P\- ■■ Pn 

On the other hand, since (^(e"*^) = ^(6*^), we see immediately that 

cos(0i + ■ ■ ■ + 0„) sin(</)i + ■ ■ ■ + 0„) 
gc\Q) = , and gs{6) = . (3.16) 

Pi--- Pn Pi--- Pn 



Therefore, using Theorem 12. 7[ with A = 1, we obtain Ii{G) — Ji{G) = ^G{0), that is 



oo 



cos(0i H h 0n) - 6'sin(0i H \- (pn) d9 _ % 

Pi---Pn l + 9^~ 2a-- ^ • ^ 

The evaluation of the integral fl3.17p . when a = 1 and tk = kr for some < r < l/ra, is 
an unsolved problem proposed by Narayana Aiyar in the begining of the twentieth century 
P], while the generalization, corresponding to a > and tk = kr for some < r < a/n, 
is a problem proposed by M. D. Hirchhorn [6j. 



Note that Theorem 12.71 yields the following more precise results, valid for A > : 

COS(01 + ■ ■ ■ + (pn) dO TT 



pi- - ■ pn A2 + 6^2 2A a - tkC 



11^-7:7^. (3-18) 



and 



sin (01 H h0n) 0d9 



Pl---Pn A2 + 02 




(3.19) 



where ti, . . . , t„ are real numbers from the interval (0, a) and the (pkS and p^'s are defined 

by dsinD. 
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Example 3.4. Bernoulli Polynomials and Polylogarithms. 

In this example we consider the sequence (-Bm)m>o of Bernoulli polynomials. They 
can be defined via the generating function 



= J2bUx) — , (3.20) 
— 1 ^-^ ml 

or recursively by -Bq = 1 and 

Vn > 0, 5; = Br,-i, and [ Bn{t) dt = 0. (3.21) 

Jo 

In particular, 

i?i(X)=X-l B2iX)=X'-X + ^, B,{X)=X{X-^){X-1). (3.22) 

Now, denote by B^ the 27r-periodic function defined by 

Bm{x) = Bm ({2"}) ' where {u} is the fractional part of u. (3.23) 

For m > 1, the Fourier series expansion of -B^, is well-known and easy to find (using 
the recursive definition fl3.2ip . see, for example [1, Chapter 23].) We have 



^ ' n=l 
n=l 



Applying Corollarv 12.61 we conclude that, for m > 1 and A > 0, we have 

r^—B fxidx-^^in^^^v^ 

/ p^^2^2'n-lWaX 2(27r)2— 2 Z^^2n^-l 

n=l 

/-°° 1 ~ _ (-l)"(2m)! >^ e- 

X p^^2r«(a;)rfx- 2;^(27r)2™-i ^ n 



^2m 

n=l 

The change of variables x -v- 2ttx and A <r- 27rA yields the following result, for m > 1 and 
A > : 

- xB,^_,{{x}) (2m -1)! 

A2 + x2 ^^-(-1) 2(271)2^-2 Li2™-i(e ) (3.26) 

(2m)! Li2^(e-^-^) 
^ A2 + x2 "^"^-^ 2(27r)2™-i A ^^'^^^ 

where the function Li^ is the polylogarithm of order k. The function Li^ is defined on 
open unit disk by the series /n^. (For an extensive account of the polylogarithms 

see El.) 
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In particular, since Bi{X) = X — i, we see that A > 0, we have 

x({x}-l/2)^^ 1,^,^ _2.A 



p + dx = -\nil-e-'-'). (3.28) 
Also, from the expressions of B2 and B3 we conclude that, for A > 0, we have 



00 



{x}(l - {x}) ^ , t: /„-27rAN 



dx = — + —- Li2(e-^'^''). (3.29) 







A2 + x2 12A 27rA 



x{x}{{x}-l)i{x}~l/2) 3 ,_2.A 







A2 + a;2 47r2 



rfx = — Li3(e~'"^). (3.30) 



Adding one forth of (13.281) to (I3.30p we see that, for A > 0, we have 

£<M^d. = - e--) + j|jLi3(e--). (3.31) 







Further Examples 3.5. 

We will end our discussion by citing some results that can be proved using the 
methods of this paper, without presenting the details : 

• For positive real numbers fi and A, we have 

xarctan(usinx) , vr , ae'^ + Wl + u? — 1 

dx = —m . 3.32 

x' + X^ 2 ^e^-^/iTJi^+i ^ ^ 

For positive real numbers ^ and A, we have 



arctan(/icosx) vr + —1 

— ax = — arctan — r . (3.33) 

x^ + X yue-^ 

For real numbers /i and A such that fi G (—1, 1) and A > 0, we have 



In l±ii^ dx = ^ln (3.34) 



a;2 + A2 1 - /i cos X A + a/1 - - 1 

For real numbers fi and A such that /i G (—1, 1) and A > 0, we have 



^ In ^ ^" — = 27r arctan ^ . (3.35) 

+ A^ 1 — /isinx jje^ 

For 2; G C and a positive A, we have 

poo iz COST TT I °° \ 

^^2^^ = X Uo(^) + 2^^"J4z)e-"M , (3.36) 

where Jn is the well-known Bessel function of the first kind of order n. 
For positive A, we have 

""^ ln(tan2 x) , vr , , , , , / „x 

, -\ — ^ dx = - In tanh A . 3.37 

x"' + A^ A 
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